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Abstract 

We propose a thermodynamic form of the equation of motion for perfect fluids 
of grade n which generaUzes the one given by J. Serrin in the case of perfectly 
compressible fluids ([1], p. 171). First integrals and circulation theorems are deduced 
and a classification of the flows is given. 
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1 Introduction 

In continuum mechanics, first gradient media cannot give a model for fluids 
with strong variations of density. Material surfaces need their own character- 
istic behavior and properties of energy |i2j. D.J. Korteweg has been the first to 
point out the convenience of fluids of grade upper than one and not to model 
interfacial layers by means of discontinuity surfaces in liquid-vapor interfaces 
[3] . Recently, this approach was used in the case of dynamic changes of phases 

S- 

Till now, the thermodynamics of fluids was neglected and it was not possible 
to model flows with strong variations of temperature such as those associated 
with combustion phenomena or non isothermal interfaces. In fact, it is difficult 
to take the gradients of temperature into account (it is not possible to consider 
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a virtual displacement of temperature) but we can use the specific entropy 
through internal energy density. 



To this aim, we usefully describe conservative flows of a compressible fluid by 
the means of its internal energy depending on both entropy and density. This 
will only be a limit mathematical model, however the study of mathematical 
structure of the equations of motion is fundamentally necessary. So, to improve 
the model for strong variations of entropy and density in interfacial layers we 
consider an internal energy function of the two quantities and their spatial 
gradients up to a r?, — 1 convenient order: the fluid will be said of grade n [S] . 

The paper aims to show that the equations of motion for perfect fluids of any 
grade can be written in an universal thermodynamic form structurally similar 
to the one given by J. Serrin in the case of conservative perfect fluids [1]. 
When the thermodynamic form is applied to the second gradient fluids [6f7] . 
it leads to three results: first integrals associated with circulation theorems 
(such as Kelvin theorems), potential equations representing the motion of the 
fluid PI5] and a classification of the flows similar to the one of compressible 
perfect fluids p^fTT] . 




2 Fluids of grade n. 

Perfect fluids of grade n {n is any integer) are continuous media with an 
internal energy per unit mass e which is a function of the specific entropy s 
and the density p in the form : 



where (grad)^, p & {1, ...,n — 1}, denotes the successive gradients in the space 
Dt occupied by the fluid at present time. We easily may consider the case 
of an inhomogeneous fluid but, for the sake of simplicity we will not do it. 
So, the material is supposed to have infinitely short memory and the motion 
history until an arbitrarily chosen past does not affect the determination of 
the stresses at present time. 



3 Equation of motion written in thermodynamic form. 

The virtual works principle (or the virtual powers principle) is a convenient 
way to find the equation of motion. For conservative motions, it writes as the 
Hamilton principle 0. 

A particle is identified in Lagrange representation by the position X(Xi, X2, X3) 
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occupied in the reference space Dq. At time t, its position is given in Dt by 
the Eulerian representation x(xi, 0:2, Xs). 

The variations of particles motion are deduced from the function family : 

X = V(x,t;a) (1) 

where a denotes the parameter defined in the vicinity of associated with a 
family of virtual motions of the fluid. The real motion corresponds with a = 

m- 

Virtual displacements associated with any variation of the real motion can be 
written in the form : 

SX = ^(x,t;a) 

da „ 
0=0 



This variation is dual and mathematically equivalent to Serrin's one p. 
145, [8]). Let L be the Lagrangian of the fluid of grade n : 

L = p Qv*V - e - 

where V denotes the velocity of particles, Q the potential of mass forces defined 
on Df and * the transposition in Dt. Between times ti and t2, the Hamilton 
action writes [T|8] : 

a = I I L dv dt. 



Iti JDt 

where dv denotes the volume element. 

The density satisfies the conservation of mass : 

pdetF = po(X) (2) 

where Pq is defined on Dq and F is the gradient of deformation. The motion is 
supposed to be conservative, then the specific entropy is constant along each 
trajectory : 

s = So(X). (3) 
Classical calculus of variations yields the variation of Hamilton action : 

From 

Sa = a'{a)\a=o, 

we deduce. 
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5a= r [ \{-- pe')6p + pViSVi- ps'Js (4) 
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The definition of dual virtual motions yields : 



dxi dx2 dxs dt. 



5(grad)^p = {gradySp and 5(grad)^s = {giadySs. 

By using the Stokes formula, let us integrate by parts. Virtual displacements 
are supposed to be null in the vicinity of the edge of Dt and integrated terms 
are null on the edge. We deduce : 



Sa — 



t2 
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Sp (5) 
Ss — pViSVi > dxidx2dxs dt. 



P< + E(-ir(P^-,n...Jm-V 



With divp denoting the divergence operator iterated p times on the edge of 
Dt, we obtain Rel. (5) in tensorial form : 



Sa 



t2 
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de 
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Ss-py*Sy} dvdt. 



By taking (2) into account, we obtain : 

<5p = pdivo5X+^||5X 

where divo denotes the divergence operator relatively to Lagrange variables in 
Do. 



We also get : 



The definition of velocity implies : 



ax, ax, , 

^(x,()V+^(x,i) = 0, 
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therefore 



Let us consider 



V + —6Y + — - = 0. 
ax ox ot 

SV = -FSX, 



where ' denotes the material derivative. Denoting 



n-l 



de 



^ = 4 + -E(-irdiVpfp 

P p=i \ 



de 
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h = e + - and m 
P 



-Y*V -h-Q, 
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then, Rel. (5) yields : 



6a - 



t2 

ti JDt 

t2 



m6p- pO 5s + p{\*Y)5X 



dv dt 



Po 



(V*F) — ^gradpS — gradgm 



SX dvn dt 



Hi JDo 

where gradg denotes the gradient operator in Dq. 

The principle for any displacement 6X null on the edge of Dq, 5a = Q implies: 

V*F = 6 gradp sq + gradg m. (6) 
(9V 

Noting that iV* + V*— )F = V*F, we get : 

ax 

F = 6'grads-grad(/i + n). (7) 
Taking Rel. (3) into account, we obtain : 

s = 0. 

Relation (7) is the generalization of Rel. (29.8) in [1]. This is a thermodynamic 
form of the equation of motion of perfect fluids of grad n. 

Obviously, term p has the same dimension as pression, 9 has the same di- 
mension as temperature and h has the same dimension as specific enthalpy. It 
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seems natural to call them pression, temperature and enthalpy of the fluid of 
grad n, respectively. 



4 Conservative properties of perfect fluids of grad n. 

Relation (7) leads to the same conclusions as those obtained in [T], [7f5|9|10|ll] 
but for fluids of grad n. Let us remind the most important results. 

With J denoting the circulation of velocity vector along a closed fluid curve 
C convected by the flow, 

at Jc 

The Kelvin theorems are deduced: the circulation of the velocity vector along 
a closed, isentropic (or isothermal) fluid curve is constant. 

For any motion of fluids of grad n, we can write the velocity fleld in the form: 

V = grad(/? + 'f/'grads + rgradx, (8) 

the scalar potentials (p, tjj, s, r and x verifying : 

^ = hy*Y -h-n, f = 0, = ^, X = 0, s = 0. (9) 

Eqations (8) and (9) induce the same classiflcation as for conservative flows 
of compressible perfect fluids PI5] : 

Oligotropic motions. - They are motions for which surfaces of equal entropy 
are vortex surfaces. The circulation of the velocity vector along a closed, isen- 
tropic fluid curve is null. Equation (8) of the motion yields : 

V = grad Lp + il) grad s. 



Homentropic motions. - In the whole fluid s is constant and Eq. (8) yields : 

V = grad + r grad x ■ 
The Cauchy theorem can be easily written : 

d rot _ aVrot V 
dt p dii. p 
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Denoting 1-L = |V*V + h + Q, Eq. (7) yields the Crocco-Vazsonyi equation 
generalized to stationary motions of perfect fluids of grad n : 

rot V X V = ^ grad s — grad 7/. 

The laws of conservation expressed by the Kelvin theorems correspond to the 
group of permutations of particles of equal entropy. 

This group keeps the equations of motion invariant. It is associated to an 
expression of Noether's theorem as in [12]. So, it is natural to conjecture such 
results for a general fluid whose internal energy is a functional of the density 
and the entropy. 
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